TROPICAL SPECTRAL CURVES AND INTEGRABLE CELLULAR 

AUTOMATA 



REI INOUE AND TOMOYUKI TAKENAWA 

Abstract. We propose a method to study the integrable cellular automata with periodic 
boundary conditions, via the tropical spectral curve and its Jacobian. We introduce the 
tropical version of eigenvector map from the isolevel set to a divisor class on the tropical 
hyperelliptic curve. We also provide some conjectures related to the divisor class and 
the Jacobian. Finally, we apply our method to the periodic box and ball system and 
clarify the algebro-geometrical meaning of the real torus introduced for its initial value 
problem. 



1. Introduction 

1.1. Background and overview. The box and ball system (BBS) [Hj and the ultra- 
discrete Toda lattice |13] are typical examples of integrable cellular automata on one- 
dimensional lattice. The key to construct these systems from known soliton equations is 
a limiting procedure called ultra- discretization [15]. These automata are also well-defined 
on a periodic lattice, and those are what we study in this paper. 

In [S], Kimijima and Tokihiro attempted to solve the initial value problem of the ultra- 
discrete periodic Toda lattice (UD-pToda). Their method consists of three steps: (1) send 
initial data of the UD-pToda to the discrete Toda lattice via inverse ultra-discretization, 
(2) solve the initial value problem for the discrete Toda lattice and (3) take the ultra- 
discrete limit. However, due to technical difficulties, this method has been completed only 
in the case of genus 1. Thereafter the initial value problem of the pBBS is solved by a 
combinatoric way and by Bethe ansatz using Kerov-Kirillov-Reshetikhin bijection [B]. 

In this paper we propose a method to study the isolevel set of the UD-pToda and the 
pBBS via the tropical spectral curve and its Jacobian intending to solve the initial 
value problem. We introduce the tropical version of eigenvector map from the isolevel 
set to a divisor class on the tropical hyperelliptic curve (Propositions 13.71 13.101 etc) . We 
provide some conjectures (Conjectures 12.31 and 13. 4|) related to the divisor class and the 
Jacobian, and also present concrete computation in the case of genus g < 3. Finally, by 
(j4.5|) we unveil the algebro-geometrical meaning of the real torus introduced in [B], on 
which the time evolution of the pBBS is linearized. 



Key words and phrases, tropical geometry, integrable dynamical system, spectral curve, eigenvector 
map, Toda lattice. 
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Tropical geometry is being established recently by many authors (see [HI [12] and ref- 
erences therein for basic literature). It is defined over tropical semifield T = MU{oo} 
equipped with the min-plus operation: "x + y" = min{x,y}, "xy" = x + y. In [9], the 
Jacobian of a tropical curve has been introduced by means of the corresponding met- 
ric graph. Our approach might be a nice application of tropical geometry to integrable 
systems and one may confirm properness of the definition in [S] . 

1.2. Tropical curve and UD-pToda. We review on how tropical geometry appears 
in studying the UD-pToda lattice. Fix g G Z>o. The {g + l)-periodic Toda lattice of 
discrete time t S Z 0] is given by the difference equations on the phase space 11 = {u* = 
(/^•••,/*+i,F/,.--,I/Vi)UeZ}c^C2(5+i): 



(1.1) 



t+i 



ft y-t 



rt+l 

i 

where we assume the periodicity Ij^g^i = and VIj^^j^^ = V^. For each n* G It, the Lax 
matrix is written as 



(1.2) 



1 



(-i)4\ 



"a 



J 



Ij+i + Vi, b\ = ijV^ and y G C is a spectral parameter. The evolution ()1.1|) 



where a* = + 
preserves det(xI+L*(y)). When we fix a polynomial f{x,y) G C[x,y] as 



+ ClX + Co) + c_i, 



(1.3) f{x,y)=y^ + y(x^+' 
the isolevel set IXc for (II. ip is 

Uc = {u'eU\ ydet{xI+L\y)) = f{x,y)}. 

Let 7c be the algebraic curve given by f{x,y) = 0. For generic Cj, jc is the hyperelliptic 
curve of genus g. Since the Lax matrix ()1.2p is same as that for the original periodic Toda 
lattice (of continuous time) [l], 11c is isomorphic to the affine part of the Jacobi variety 
Jac(7c) of 7c, and the time evolution (II. Ih is linearized on Jac(7c) [5]. 

The ultra- discrete limit of (jl.ip is the UD-pToda [11] given by the piecewise- linear map 

{t £ Z and i,j G {1,2, . . . , g + 1}), where 

(1.4) Q*+i = min[W,*, Qj - Xj], W^+' = Q^, + - Q'+\ 



0, 



T 1 1 1 — 

Al A2 • • • Xg-l \g 



Figure 1. Tropical hyperelliptic curve 

with Xj = minfc=o,--- ,g [Sf=i(^/-/ ~ Qi-z)] • other hand, in this hmit 7c is reduced 

to the tropical curve Tc C given by the polygonal lines of the convex in M^: 

(1.5) {iX,Y,mm[2Y,{g + l)X + Y,gX + Y + Cg,--- , X + Y + Ci,Y + Co,C^i])} 

For generic Ci (see (j2.ip ) , Tc is smooth and depicted as Fig. [1] where we fix = and 
set Aj = Cg-i — Cg-iJ^i for i = 1, • ■ ■ ,g. Note that all edges of F have fractional slopes. 

We explicitly construct a tropical version of the eigenvector map from the isolevel set 
of the UD-pToda to the divisor class on Fc, and show that the isolevel set is isomorphic 
to the tropical Jacobi variety of Fc*. 

1.3. Ultra-discrete limit and min-plus algebra. We briefly introduce the notion of 
the ultra-discrete limit (UD-limit) and relate it to the min-plus algebra on the tropical 
semifield T = ]RU{oo}. 

We define a map Log^ : IR>o — > M with an infinitesimal parameter e > by 

(1.6) Log;, : X I— > —elogx. 

For rr > 0, we define X G T by rr = e~~. Then the limit e — > of Log£(x) converges 

to X. The procedure limE^oLogg with the scale transformation a.s x = e~~ is called the 

ultra-discrete limit. 

We summarize this procedure in more general setting: 
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Proposition 1.1. For A,B,C G M and ka, kt, kc > 0, set 

a = kae~~, b = kf)e~~ , c = kce~~ 
and take the limit e ^ of the image Log^ of the equations 

(i) a + b = c, (ii) ab = c, (iii) a — b = c. 

Then 

(i) min[^, B] = C, (ii) A + B = C 

and 

( A = C (ii A< B, or A = B and ka > kb) 

(iii) { 

I contradiction (otherwise) 

hold. 

1.4. Content. In §2, we define the metric graph Tc for the tropical hyperelhptic curve 
Tc and define its Jacobian JiVc)- By using a tropical version of the Abel-Jacobi map, we 
propose a divisor class which isomorphic to JiJ^c) at Conjecture 12.31 This is justified for 
g( < 3. In §3, we study the isolevel set of the UD-pToda. We construct the eigenvector map 
from the isolevel set to the divisor class on the tropical curve. It is shown that the general 
level set is isomorphic to J^Tc)- In §4, we clarify the correspondence of the UD-pToda 
with the pBBS by refining that in [3]. In conclusion we interpret the isolevel set of the 
pBBS introduced in [5] in terms of tropical geometry. 

Acknowledgement. R. I. thanks Atsuo Kuniba for discussion. She also thanks Keigo 
Hijii for quick help in using LaTeX. T. T. appreciates the assistance from the Japan Society 
for the Promotion of Science. R. I is supported by Grant-in- Aid for Young Scientists (B) 
(19740231). 

2. Tropical hyperelliptic curve and Jacobian 

2.1. Tropical hyperelliptic curve. Fix g G Z>o and C = (C_i,Co,-- - ,Cg) G M^^^ 
Let Tc C M? be the affine tropical curve given by the polygonal lines of the convex in 
(|1.5|) . We assume a generic condition for C: 

(2.1) C_i > 2Co, Ci + C,+2 > 2Q+1 (i = 0, • ■ ■ ,9 - 2), Cg.i > 2Cg. 

For simplicity, we fix Cg = in the following. Define A = (Ai, • • • , \g) and pi, • • ■ ,Pg by 

9 

(2.2) \i = Cg-i- Cg-i+i, Pi = C-i - 2y^min[Ai, Aj]. 

i=i 

Under the condition (|2.ip one sees < Ai < A2 < • • • < Ag and 2 J2i=i -^j < C_i. 

By referring [U Definition 2.18], we introduce a notion of smoothness of tropical curves: 
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Definition 2.1. The tropical curve S is smooth if the fohowing conditions are 

satisfied: 

(a) all edges in T, have fractional slopes. 

(b) All vertex G S is 3-valent. 

(c) For each 3-valent vertex v, let 61,62,63 be the oriented edges outgoing from v. Then 
the primitive tangent vectors of satisfy Ylt=i^k = 0, and j^fc /\ Cj\ = 1 for k / j, 
fc,iG {1,2,3}. 

We see that Tc is smooth. In particular, it is a tropical hyperelliptic curve whose 
genus is dim Hi {Tc = g (see Fig. [1]) . We are to consider the maximal compact subset 
Tc = Tc \ {infinite edges} of Tc- For simplicity we write F for Tc- 

2.2. Metric on F. Following [9l §3.3], we equip F with the structure of a metric graph. 
Let £(F) be the set of edges in F, and define the weight w : £(F) M>o by 

II ?e II 

where is the primitive tangent vector of e G £(F), and || || denotes any norm in M."^. 
With this weight the tropical curve F becomes a metric graph. 

The metric on F defines a symmetric bilinear form Q on the space of paths in F as 
follows: for a non-self-intersecting path 7, set (5(7,7) := length^(7), and extending it 
to any pairs of paths bilinearly. In Fig. [2] we show the weight for each edge in F and 
the basis (i = 1, • • • , of 7ri(F). For example, we have Q{ai,ai) = C_i + pi + 2Ai, 
Q{ai,a2) = —pi, and Q{ai,ai) = for i > 2. 

2.3. Tropical Jacobian. Let fl{T) be the space of global 1-forms on F, and r2(F)* be the 
dual space of 0(F). Then both 0,(T) and r2(F)* are g dimensional and r2(F)* is isomorphic 
to Hi{T,R). 

Definition 2.2. [9l §6.1] The tropical Jacobian of F is a 5 dimensional real torus defined 
as 

J(F) = n{Ty/Hi (F, Z) ~ R3 /K Z3 ~ /a . 
Here K,A£ Mg (M) are given by 

Kij =Q{ai,aj), 

i j 

Kj =QC^ak,'^ai) = C-i +Pi5ij + 2m.m[Xi,Xj]. 

k=l 1=1 

Since Q is nondegenerate, K and A are symmetric and positive definite. In particular, 
we say that J(F) is principally polarized. 
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\o — Ai 



V2 





Figure 2. Tc as a metric graph 



Let Div|^g(r) be a set of effective divisors of degree g on T. We fix Pq ^ T and define a 



map T] : Div^g(r) ^ J(r); 
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(2.3) Pi + --- + Pg^ Y^^Qili, ai), • ■ ■ , g(7„ aj), 



where 7i is the path from Pq to Pj on F. Define Ojj = ajnaj\{the end-points of Haj} C 
F, and ©^(F) to be a subset of Div^g(F): 



'D5(F) 



A + • • • + n 



Pj G Oj for all i, and 

there exists at most one point on aij for all i ^ j 



Conjecture 2.3. A reduced map r/|D9(r) bijective: 

7?|a,«(r) : 2)^(r) ^ J(F). 

In the case of g = 1, this conjecture is obviously true since D^(T) = F ~ >/(F) by 
definition. In the following we show that this conjecture is true for g = 2 and 3. 



Proof. We define a map is ■ T 



P ^ tsiP) = {Qil,Oii))i<i<g where 5 G F and 7 is 



an appropriate path from S to P. For Pi + ■ • • + E DiVgg(F), we see 7?(Pi + • • ■ + Pg) ~ 
ELi ^Po(^'*) in J(F). 

g = 2 case: We set Pq = (Ai,2Ai) which is the end-point of 012- In the left figure of 
Fig. Owe illustrate the locus of tP(, (Pj) where Pj starts from Pq and moves along Oj for 
i = 1,2 respectively. We set O = (0,0), Ai = (C_i + pi + 2Ai, -pi) and A2 = {-pi,2pi). 
The parallelogram F of dash lines is the fundamental domain of J(T). We calculate the 




Figure 3. Imege of ij in g = 2 



image of the map ©^(r) given by Pi + P2 ^ '-Po(-Pi) + ^Po{P2) + M0-, and obtain 

the parahelohexagon V composed of three non-overlapped parallelograms as shown in the 
right figure of Fig. [3l It is easy to see that V is isomorphic to F in JiT). 

g = 3 case: We set Si = (Ai,3Ai) and ^2 = (A2, Ai + 2A2), the end-points of au and 
023 respectively. We calculate the image of the map D^(r) — s- M.^ given by Pi + P2 + P3 1-^ 
i-SiiPi) + i-Si{P2) + i-S2iP3)j and obtain non-overlapped 12 parallelopipeds. After shifting 
some parallelopipeds along the lattice Kl?, we obtain the parallelo dodecahedron V in 
Fig.H We set = (0,0,0), ^1 = (C_i + pi + 2Ai, -pi, 0), ^2 = (pi, -pi - P2 - 2(A2 - 

A2),P2), ^3 = (0, -P2, 2p2), P = ^1 + ^2 + ^3 and Bi = Aj + Ak for {i,j, k} = {1, 2, 3}. 

> > > 

The parallelopiped F spanned by OAi,OA2 and OA^ is the fundamental domain of J{T). 
We draw V in black, and F in blue. 

One sees that V coincides with F in J(r) as follows: The polygon V \ F is composed 
of three parts each of which contains the face *' in Z2-23-plane, the face * in zizs-plane or 
the face o' in ziZ2-plane. We translate the part with the face *' (resp. o') by OAi (resp. 
OA2, 0^43) and attach it on the face * (resp. o). □ 

Remark 2.4. After this paper was submitted, we proved Conjecture 12.31 for general g in 
another way, by applying the notion of rational functions on F [31 [9] . 

For § 4.2, we introduce another torus J'(F): 

(2.4) J'{T) =R9 /AZ^, 

where Aij = Aij — C_i. 

Lemma 2.5. (i) det A = detK = {g+ 1) detA = {g + l)pi ■ --pg-iC-i. 

(ii) Let u-p he a shift operator, up : — > M^; (zj)j=i,... ,g 1— > (zj -|- C_i)j=i^... ^g. Then 
J'(F) ~ J(F)/{P - ur{P) I P G J(F)}. 

The proof is elementary and left for readers. 
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3. ISOLEVEL SET OF ULTRA-DISCRETE PERIODIC TODA 

3.1. Periodic Toda lattice. We review the known results on the {g + l)-periodic Toda 
lattice (ll.ip . We define a matrix M^{y) besides the Lax matrix L^{y) ()1.2p on the phase 
space U: 



/4 1 



M*(y) 



1 

-'3 ^ 



Proposition 3.1. [3] (i) The system (jl.ip is equivalent to the Lax 



orm 



L'+\y)M\y) = M\y)L\y) 



(ii) The system (jl.ip preserves the characteristic polynomial of L*{y), det(x I +L*(y)). 
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Proof, (i) Set R\y) as 



(-1> 



vi 1 



V 



J 



The system ([TI]) is equivalent to i?*+^(y)M*+i(y) = M\y)R\y). By the fact 
R\y)M\y), we have 



L*+^(y)M*(2/) = i?*+^(y)M*+^(y)M*(y) = M*(y)i?*(y)M*(y) = L\y)M\y). 

(ii) From the Lax form we obtain det(x I = det(xI+M*(?/)L*(y)(M*(?/))~i) = 
det(xI+L*(y)). □ 

We define the (complex) spectral curve 7c given by 

f{x,y) =ydet{Ix + L\y)) 

=y2 + y(x^+^ + Cgx^ + • • • + Co) + c_i = 0. 

Concretely, Cj is given by (for simplicity, we write /* = Ii,V^ = Vi and so on) 

^9= J2 E 



(3.1 



(3.2) 

5+1 5+1 

j=i j=i 

9+1 
i=l 

For generic Cj, 7c is a hyperelliptic curve. Since (II. ip is invariant under (/j, Vi)i<j<g+i i— > 
(IjCg, ViCg)i<i<g+i, we can set Cg = 1 without loss of generality. 

Proposition 3.2. [5] Under the condition Y[k=\ W^k=i ^1' system (jl.ip is equiv- 
alent to the system: 



Tt+l 



1 _ n^^+^ % 



(3.3) 



1 _L - 



t+i 



ji+1 • 



T 

i — k 



3.2. Ultra-discrete Toda lattice. Suppose 

Vl > 0, Z* > 0, 



(3.4) 



ni'.'<n'^- 

i=l i=l 



Qi 



In the UD-limit Wm^^Q Log^ with the scale transformation li = e ^ ,Vi = e e, the sys- 
tem (j3.3p becomes the UD-pToda lattice (|1.4p . Simultaneously, the limit of the conserved 
quantities Cj = e ^ become 

Co = min[ min Qi, min VFjl, 

C„_i = min[ min (Qi + Q,), min (Wi + Wj), min (Qi + W,)], 



Co = mm[^Qi,^Wi], 

i=l i=l 

9+1 

C_i = ^(Q,+T^i), 

i=l 

which are preserved under (II. 4|) by construction. Prom the assumption (13.40 . we have 

9+1 9+1 
j=l i=l 

We can set Cg = without loss of generality corresponding to Cg = 1. 

Next, we reconstruct the tropical curve Tq by the UD-limit of the real part of the 
spectral curve 7c. We write 71R for the real part of 7 = jc- Then the image of the map 
Log^ : — > M^; (x, y) 1— > (log log \y\) of 7k is the rim of the amoeba of 7. 

In taking the UD-limit of the equation (j3.1|) with the scale transformation Cj = 
e~~^ , \x\ = e~T and \y \ = e~~, we have the following: 

(i) X > 0, y > leads to a contradiction. 

(ii) X < 0,y > 0. We have 

min[2y, C7_i, {g + l)X + Y,{g - 1)X + Y + Cg.i, . . . ,Y + Co] 
= inm[gX + Y + Cg,{g - 2)X + Y + Cg_2, ...,X + Y + Ci] 
min[2y, C.i,gX + Y + Cg, {g-2)X + Y + Cg_2, . . . ,Y + Co] 
= mm[{g + l)X + Y,{g - 1)X + Y + Cg.i, ■ ■ ■ ,X + Y + Ci] 

(iii) X < 0, y < 0. We have 



r2 : < 



mm[2Y,C.i,gX + Y + Cg,ig - 2)X + Y + Cg_2, ...,X + Y + Ci] 
= mm[{g + 1)X + Y,{g-l)X + Y + Cg.i, . . . ,Y + Cq] 
mm[2Y, C_i ,(g + l)X + Y,{g - 1)X + Y + Cg.i, . . . ,X + Y + Ci] 
= mm[gX + Y + Cg,ig-2)X + Y + Cg.2, . . . ,Y + Co] 



[g : odd) 
{g : even) 

{9 ■■ odd) 
{g : even) 
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(iv) X > 0,y < 0. We have 

r4 : min[2y, = mm[{g + 1)X + Y,gX + Y + Cg, 
Then we obtain the following. 
Proposition 3.3. For generic Ci 's which satisfy 12. 

f c = r2 U Ts = U Ta U r4 

hold. 

Fig. shows 7m, T3 and r4 in the case oi g = 2. 



o 



-1 




-3 -2-10 1 
(C_i,Co,Ci,C2) = (20,7,2,0) 
e = 4 



C_i 



Ai A2 



C_i 



Figure 5. Real and tropical curves 
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3.3. Eigenvector map. Let T be the phase space of the ultra-discrete {g + l)-periodic 
Toda lattice, and C be the moduli space of the compact tropical curves Tq- 

9+1 9+1 

7={{Qir-- ,Qg+i,Wi,--- ,Wg+i)eR''3+'' I J^Qi<5^Ty.}, 

i=l 1=1 

e = {c = • • • , Cg.i,Cg)} ~ . 

We define a map $ : 1 ^ 6 by i^, and set 7c = '^"^{C) for C G C. 

Conjecture 3.4. For a generic C = (C_i,--- ,Cg-i,Cg = 0) € S which satisfies (12. ip . 
following are satisfied: 

(i) Tc~ J(rc). 

(ii) Suppose C G Z^"*"^, and /ei (Tc')z Q^'^f^ <^z(rc) ^/'■e seis o/ lattice points in 7c and in 
J{Tc) respectively. Then the isomorphism of (i) induces the bijection between (7c)z and 
Ji^i^c)- In particular, we have \{7c)z\ = det A. 

Remark 3.5. This conjecture claims that we need only a compact part Tc of Tc to 
construct the isolevel set 7c- 

In the rest of this section, we construct the isomorphism vr : 7c — ^ J(^c) in the 
case of g = 1,2 and 3, by applying the technique of eigenvector map, which is essentially 
the same with Sklyanin's separation of variable in our case (for example see [21 [TO]). The 
isomorphism vr is a composition of isomorphisms: 

7c ^ V^iTc) ^ J(rc), 

where is called the eigenvector map (or separation of variables) and r] is the Abel-Jacobi 
map 

Remark 3.6. By concrete computation we also conjecture the following. Define a trans- 
lation operator v as 

V : J{Tc) J(Xc); z + (Ai, A2 - Ai, . . . , Ag - \g-i). 

Then the following diagram is commutative: 

Tc ^ J{Tc) 
It i V 

Tc ^ J(rc) 

i. e. the flow of the UD-pToda is linearized on the tropical Jacobian. It is easy to check 
it in the case of = 1. 

First we discuss the discrete case. Let us consider the eigenvector (p of the Lax matrix 
L^{y)- Then (j) is given by 

(k = *(/l; /2, • • • , /g, -fg+l), 
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where fi{i = l,2,...,g) is 



fi = det 



1 2 

hl+X li2 

hi I22 + X 



h,g+l 



and 



\ ^3,1 ^3,2 ••• ^ 

fhl+X li2 

hi h2 + a; 



3,3+1 



g 

h,g 
h,g 

h,g + ^ / 



fg+i = det 



3,2 



h,g 
h,g 

^3,3 + ^ 



where = {L*{y))ij. The equation fg^i{x) = has the solution xi,a;2, • • • ,Xg, each of 
which defines two points on {xi,yi), {xi,y'-), where one of them (we assume that is 
{xi,yi)) leads /,• = for all j. We choose {{xi,yi) \ i = l,2,...,g} or | i = 

1,2, ...,(/} as a representative of Pic^(7c). In the discrete case, this map induces an 
injection lie ^ Pic^(7c)) and the evolution equation (II. Ih is linearized on the Jacobi 
variety of 7^, Jac(7c) ~ Pic5(7c) (Cf. [DEI [TO]). 



3.4. The case of g = I. The Lax matrix is 



ai 



y 



■ y a2 



and the conserved quantities are 

c_i = 6162, Co = aia2 - bi — 62, ci = ai + 02- 

When /2 = ai + X = 0, (13. ip becomes 

f{x,y) = {y-bi){y-b2) = 0. 

Thus we define the map Uc 7c by n* i-^^ (xi = — ai, yi = 61). 
In the ultra-discrete limit, the map tjj : 7c ^ is given by 

{Qi,Q2, Wi, W2) ^ {Xi = min[Q2, VFi], >1 = Qi + Wi) € r2 

= Qi + (52 and Ci = min[Qi, Q2, W^i, W2] 



where C_i = Qi + Q2 + ^^i + W2, C, 
see the following: 



0. We 



Proposition 3.7. The map ip is bijective. In particular, 7c ~ J{Tc)- 
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Proof. By construction it is obvious that the image of ip is included in Tq- Inversely, 
solving 



ai = -X, bi = y 

Co + 61 + 62 , c_i 

02 = , 62 = 

ai y 

for Ij, Vj, we have the solutions (/j, Vj) and (I^', T^') = 1, 2): 

, co + 2y , x(2c_i + Coy) 

^1 + ^1 = , i2 + -'2 = : : — 2 

X c-i + Coy + 

^1 A yi2 

where we assume /j > I'^. Only (/i,Vj) satisfies the assumption (I3.4p . By the UD-limit, 
we have the inverse of ip as 

Qi =min[Co,y] -X 

Q2 =X + min[C7_i, Co + Y] - min[C_i, Cq + 2y] 

Wi =Y - Qi 

W2=C-i-l^-(32. 

□ 

3.5. The case oi g = 2. In this and the next subsection we denote min[ ] simply by 
[ ]. The Lax matrix is 



and the conserved quantities are 



c_i = 6162&3, Co = 010203 - 0261 - 0362 - 0163, 

Cl = 0102 + 0203 + 0301 — 61 - 62 — C2 = ai + 02 + 03. 

The UD-pToda (II3D is 

.3g^ Q't' = \Wl,Qt-Xl] 



with 



X* = [0, wti - Qli, wU + wU- Qi-i - QUI 
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and the conserved quantities ()3.5p become 

C2=[Ql,Q2,Q3,M^l,W^2,VF3] =0 

Ci =[Qi + Q2, Q2 + Q3, Qi + Qi, Wi + TVa, W2 + W^3, W^ + Wi, 

Qi + W2,Q2 + W^,Q^ + Wi] 
Co =[Qi + Q2 + Qs, W1 + W2 + Ws] 
C_i =Qi + Q2 + Q3 + W1 + W2 + W3. 

The tropical spectral curve is the set sum of 

: [2Y,C-i,2X + Y + C2,Y + Co] = [3X + Y,X + Y + Ci] 

and 

T3 : [2Y,C^i,3X + Y,X + Y + Ci] = [2X + Y + C2,Y + Co]. 
The eigenvector of the Lax matrix is given by 



/i 



h. 


1 


, /2 = 


ai + X 


h. 


1 h — 


y 




y 


1 


a2 + X 




62 


1 





ai + X 1 
62 a2 + X 



When /s = 0, ([Si]) reduces to 

f{x, y) = {y- bi{x + a2)){y - 63(2; + ai)) = 0. 

The solutions are 

xi+X2 = -ai - 02, X1X2 = aia2 - 62 

= + 02), y[ = b2,{xi + ai) for ? = 1,2. 
For the UD-limit we use another expression of yi: 

y 



(3.7) 



C-l 



bsixi + ai)' 

Under the assumption xi,X2 < 0, yi < 0, 7/2 > (for small e > 0), the UD-limit of 



Xl =[Q2,Q3,Wi,W2] 

X2 =[Q2 + Qs, Wi + W2, Qi + Wi] - Xi 



(I37n) leads: 
(3.8) 

and 



Yo 



Yf := Q1 + W1 + Xi if Xi < [Q3, W2] 

:=C^i-iQ3 + W3 + Xi) if Xi < [Q2, VFi] 

Y^- := Q1 + W1 + [Qs, W^2] if X2 > [Q3, W2] 

y| := C_i - (Q3 + W3 + [Q2, W^i]) if X2 > [Q2,Wi] 



The following lemma can be proved elementarily. 
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Lemma 3.8. (i) C2(= 0) < Xi < C2 + Ai < X2 < C2 + A2. 

(ii) Xi = [[Q2,Wi],[Q3,W2]], X2 > max[[Q2,Wi],[Q3,W2]]. 

(iii) If[Q2,Wi] = [Q3,W2], then 

(iii-l)Xi = X2 and thus = and = Y^ 

or (iii-2) Yf = Y^ and Y^ = Y^ 

hold. 



By Lemma 13.81 the correspondence between {Qi,Q2,Q3,Wi,W2,Ws) € 7c and 
{Xi,Yi) + {X2,Y2) G Div^g(rc) is uniquely expanded as a continuous map : Tc — > 
Div2g(rc) by ^ and 

Y2 = Y^^ = C-i-(Q3 + W3 + lQ2,W,]) J ' 
Y2 = Yf = Q,+Wi + lQ3,W2] J ' 



(When Xi = X2, we can exchange Yi and Y2.) 



Lemma 3.9. The image of ip is included in 'D'^{Tc), i-e. if Xi = X2, then {Xi,Yi) or 
(X2,y2) is at the end point 0/012- 

Proof. By Lemma [3.8l fii). we have [(52iW^i] = [Q3>W^2]- Without loss of generality we 
can assume Qi = 0. (i) Q2 = Qi Wi,W2 leads Ci = Q2 and C2 = 2Q2, which 
contradict to the smoothness ()2.1|) . (ii) Q2 = W2 < QsjWi leads Xi = Q2 and X2 > Q2', 
which is a contradiction, (iii) Wi = Q3 < Q21W2 leads Ci = Wi and Y^ = 2C'i. (iv) 
Wi = W2<Q2, Q3 leads Ci = Wi and Yf = 2Ci. □ 



Inversely, solving 



ai 



a2 



as 



xivi - X2y2 

yi-y2 

_ xiy2 - X2yi 

yi-y2 

coai&s + 02^1 

0102 — 62 



yi - j/2 

Xi - X2 

ymjxi - X2f 

{yi - y2Y 
c-i{yi - j/2) 
yiy2(xi - X2) 



for Ii, Vj, we have (e.g.) 



J ^j' ^ co(xi - X2) + 2{xiy2 - X2yi) 

^ XlX2(3;i - X2) 
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By the UD-limit, we have the inverse of ip if Xi < X2 
Qi=[Co + Xi,C/2]-(2Xi + X2) 

Q2 =2X1 + [C_i + f/i, + ^2 + f/2, Co + [Xi + Yi + Y2,X2 + 2[Yi , Y2]]] 

- [Yi,Y2] - [C_i + 2X1, Co + Xi + [/2, 2U2] 

Q3 =Xi + X2 + [yi, I2] + [c_i + f/i, Co + Xi + + ^2] 

- [C_i + 2C/i, Co + Xi + Yi + + Ui,2Xi + 2^ + 2^2] 
Wi = [Yi,Y2] - Xi - Qi 

W2 = Yi + Y2 + 2X1 - 2[Yi,Y2] - Q2 

= [C_i + [^1,^2] -Yi-Y2-Xi-Q^ 

with 

C/i = [Xi +Yi,X2 + Y2I U2 = [Xi + Y2.X2 + Yi]. 

By Lemma 13.91 the inverse is uniquely expanded as a continuous map to the case of 
Xi = X2. 

Now we have the fohowing. 

Proposition 3.10. The UD-eigenvector map ip : 7c T)'^{Tc) is bijective. 

3.6. The case of g = 3. In the case of 5 = 3 we present the ultra-discrete eigenvector 
map Tp : 7c D^{Tc)- However, for the reason of complexity, we will omit to present 
the inverse mapping and to prove the bijectivity. 
The solutions of /4 = and f{x, y) = are 

xi + X2 + X3 = -ai -02-02 

X1X2 + X2X3 + X3X1 = 0102 + 0203 + 0301 -62-^3 
X1X2X3 = -010203 + 01O3 + 03O2 

yi = fei((o2 + 2;i)(o3 + Xi) - 63)), y'i = ^{{ai + Xi)(o2 + Xi) - 62)) for i = 1,2, 3. 
For the UD-limit we use other expressions of yi: 

■ - b b ^' - £zi 

^oi + Xj b4{{ai + Xi){a2 + Xi) - b2))' 

The UD-limit of jSJ]) leads the UD-eigenvector map iIj : 7c ^ ©^(rc): 

Xl =[Q2,Q3,Q4,Wi,W2,W3] 

X2 =[Q2 + Q3, Qs + Q4, Q2 + Q4, Wi + W2, W2 + W3, Wi + W3, Qa + Wi, 

Qa + W2, Q2 + W3, Q3 + Wi] - Xl 
X3 =[Q2 + Q3 + Q4, W1 + Q3 + Qa, T^i + P^2 + Qa, W1 + W2 + W3] - {Xi + X2) 
Y,=Y^^ fori = 1,2, 3, 
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where 



Yl =Qi + Wi + [2Xi,Xi + [Qs, Q4, W2, W3], [Q3 + Q4, W2 + W3, Q4 + W2]] 

y.2 + W1+Q2 + W2+ [Q4, W3,X,] - [Q2, Wi, X,] 

=c„i - (Q4 + + [2X„ + [Q2, Q3, Wi,W2], [Q2 + Q3, Wi + i^2, Qs + Wi]]) 

and Sj is defined as follows. 

(i) Set A\A'',A\B\B^ as = [g2,VFi],^' = [Q3,W2],A^ = [Qa,W3], 

= [Q3 + Q4, W2 + W3, g4 + W2], 53 = [Q2 + Q3, W^i + W2, Q3 + Wi], and define si by 

si = 1 if ^1 < [A^,A% 
si = 2 if ^2 < [^3^^1]^ 

51 = 3 if ^3 < [^1,^2]^ 

If si has two or more possibilities, choose one of them. 

(ii) Define S2 and S3 so that Sj 7^ Sj{i,j = 1,2,3) by 

52 = 1 if X2 + [^2^ 

S2 = 2 if < X2 < A^ or ^3 < X2 < 

52 = 3ifX2 + [A\A^]<B^ 
and 

53 = 1 ifX3 + [A2^A3] 

S3 = 2 if X3 > max[A\A3], 
S3 = 3 if X3 + >^3_ 

(iii) If both S2 and S3 are not determined by (ii), then choose S2 and S3 arbitrarily under 
keeping Sj / Sj{i,j = 1,2,3). 

4. From the UD-pToda to the pBBS 
4.1. The structure of T^. Fix a generic C G C with Cg = 0. Define a shift operator 

(4.1) (Qi,... ,Qg+i,T^i,... ,W,+i)^ (Q2,-- - ,Qg+i,Qi,W2,--- ,Wg+i,Wi). 

Note s^"*"^ = id. We define a subspace of 7c- 

(4.2) 

7c = {(Qir-- ,Wg+i)e7c I (a) I^i > 0, and (b) Qi = or Wg+i = 0.} . 

We write Ti^^ for the set given by 

= {s\t) I r e T^}, for i G Z. 

Proposition 4.1. (i) D = for i ^ j mod 5 + 1, (ii) 1c = y T^. 

First we show 
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Lemma 4.2. If t = (Qi, • • • , Qg+i,Wi, ■■■ , Wg+i) G T^, then Qi > for 2 < i < g, and 
Wj >0 forl<j <g. 

Proof. Recall that the conserved quantity Cg_i (|3.2|) for Tc satisfies Cg_i > 0. For 
T = (Qi, • ■ ■ 1 Qg+i, Wi, • • • , Wg+i) € Tc, the following properties (bl) and (b2) hold: 
(bl) When Qi = 0, we have 

(4.3) C„-i = min[ min Qi, min Wi,Wi + W„+i\ > 0. 

Thus we obtain Qi > for 2 < i < g + 1 and Wj > for 2 < i < g. 
(b2) When Wg+i = 0, we have 

(4.4) C„_i = min[ min Qi, min Wj, Qi + Qg+i] > 0. 

Thus we obtain Qi > for 2 < i < g and Wj > for 1 < i < g. 

If we further assume r G T^, we have Wi > 0, and (bl) or (b2) is satisfied. Thus one 
obtains the claim. □ 

Proof. (Proposition 14. 1|) 

(i) Note that 

Th = {{Qu ■■■ , Qg+uWi, ■■■ , Wg+i) I (a) Wi+i > 0, and (b) Q,+i = or W = 0}. 

We check that if r = {Qi, ■ ■ ■ , Qg+i, Wi, • • • , Wg+i) G then it satisfies (a') Wj+i = 0, 

or (b') Qi+i > and Wj > 0, for i = 1, • • • , (7. For i = 1, - ■ ■ ,g — 1, (b') is satisfied due to 

Lemma |4.2[ For i = g, (b') is satisfied when Qi = and (a') is satisfied when W^+i = 0. 

9 9 

(ii) Is is trivial that Tc D |J ^c- We show 1c C |J T^. Since Cg = 0, for r G Tc we 

1=0 i=0 

assume Qi = or W^+i = without loss of the generality. When Qi = 0, (14. Sp denotes 
Q2, • • • , Qg+i, W2, • • • , Wg > and Wi + Wg+i > 0. Thus we see t e T^, when Wi = 0, 
and r G when Wi > 0. In the same way, when W^+i = it is easy to see that ()4.4p 
indicates r G T^. □ 

4.2. Periodic BBS. FixLGZ>o. The L-periodic box-ball system is a cellular automaton 
that the finite number of balls move in a periodic array of L boxes, where each box has 
one ball at most [16]. We assume that the number of balls |A| satisfies 2|A| < L. The time 
evolution of the pBBS from the time step t to t + 1 is given as follows: 

(1) Choose one ball and move it to the leftmost empty box to its right. 

(2) Choose one of unmoved balls and move it as (1), ignoring the boxes to which and 
from which the balls were moved in this time step. 

(3) Continue (2) until every ball moves once. 

This system has conserved quantities parametrized by a non-decreasing array A = 
(Ai, • • • ,Xg) G (Z>o)^ such that Yli=i — \M for some g G Z>o. In the following we write 
and 1 for "an empty box" and "an occupied box" respectively, and let Bl ~ {0,1}^^ 
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be the phase space of L-periodic BBS. We show examples of the evolution of b{t) € as 
time t goes: 

Example 4.3. The case of (i) (L, Ai) = (8, 3) and (ii) (L, Ai, As) = (7, 1, 2): 







(ii) 








+ 

L 




L 


0{l) 




1 [P{0{U))) 





00111000 





0100110 


(0,1,2,1,2,1) 


(0,1,2,1,2,1) 


1 


00000111 


1 


1010001 


(1,1,1,1,3,0) 


(1,1,1,1,3,0) 


2 


11100000 


2 


0101100 


(0,1,2,1,1,2) 


(1,2,0,1,2,1) 


3 


00011100 


3 


0010011 


(0,1,2,2,2,0) 


(1,2,0,2,0,2) 


4 


10000011 


4 


1101000 


(2,1,0,1,3,0) 


(1,0,2,3,0,1) 


5 


01110000 


5 


0010110 


(0,1,2,2,1,1) 


(2,0,1,1,2,1) 



Roughly speaking, g is the number of consecutive clusters of I's, and (Ai, • • • , Xg) cor- 
responds to the number of I's in each cluster. 

The injection from Bl to T is introduced in [5]. Fix a generic C e en which 
satisfies (f2T]l with C_i = L and Cg = 0, and set A = (Ai, • • • , Xg) ([221). Note that the 
generic condition for C corresponds to the condition: < Ai < A2 < • • • < Ag. Let 
Bi x C Bl be a set of the states whose conserved quantity is A. Then the injection 
P ■ Blx ^ ('J'c)z; b (Qi, • • • , Qg+i, Wi,--- , Wg+i) is defined as follows: 

(1) if the leftmost entry of b is 1, then set Qi = tJ(the first consecutive I's from the left), 
otherwise set Qi = 0. 

(2) Set Wi = ti(the i-th consecutive O's from the left) for i = 1, ■ ■ ■ , g + I. If 
Qi 7^ 0, set Qi = tj(the i-th consecutive I's from the left), otherwise set Qi = 
ji(the (i — l)-th consecutive I's from the left) for i = 2, • • • ,g + l. 

Proposition 4.4. /3 : B^ x (^c)^ ^c* ^ '^'^^^'^^^ is a bijection. In particular, 
we have the bijection between {7c)z/{t ~ s{t) \ r E {7c)z} CL^d B^^x, which leads to 
l(Tc)z| = {g + l)\BLA\. 

Proof. By the definition of the map (3, it is obvious (3{Bi x) C {T^)z- From Lemma 14.21 
each T gives the array {Qi, Q2, • • • , Wg, Qg+i,Wg+i) where Wi,Q2, ■ ■ ■ ,Wg > 0. 
We define a map p : {Tq)i — > B^^x as follows: when Qi = 0, set p{t) as 

0_^1_^ ■ ■ ■ l_^Oj^ 

Wl Q2 Qg + l Wg+1 

where Wg+i can be zero. When VFg+i = 0, set rj{T) as 

1_^0_^ . . . 0_^1_^ 

Ql Wl Wg Qg + l 

where one of Qi and Qg+i can be zero. In both cases, it is clear that (3 ■ p{t) = r. Thus 

p = p-\ □ 
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From Prop. STT] and [131 we can put back b{t) from the solution of the UD-pToda lattice 
with the initial state /?(6(0)) (see Example 14.31 (ii)). 

Lemma 12.51 clarifies the algebro-geometrical meaning of J'{Tc) ()2.4p which was first 
introduced in the study of the pBBS by Kuniba et al: 

Theorem 4.5. [U Theorem 3.11] Let j!^{Tc) be the set of lattice points in J'{Tc)- Then 
the bijection between Bj^ x and Jl^iTc) is induced by Kerov-Kirillov-Reshetikhin bijection. 

In the following diagram we summarize the relation among the UD-pToda, the pBBS 
and the tropical Jacobian: 

/3 

Bl,x ('J'c)z C 7c 

(4-5) ii ^ li ii 

Wc) -/.^ MTc) C J{Tc) 

Here /s and /tr are the quotient maps respectively induced by the shift operators s (14. ip 
and at Lemma [23] (ii). The isomorphism of the right two downward maps are the claim 
in Conjecture [M The diagram also indicates J'{Tc) ^ J(Tc)/{P ~ s*{P) \ P G J{Tc)]- 

References 

[1] M. Adler and P. van Moerbeke, Linearization of Hamiltonian systems, Jacobi varieties and represen- 
tation theory, Advances. Math. 38, 318-379 (1980). 

[2] M. Audin, Spinning tops - A course on Integrable Systems, (Cambridge University Press, 1996). 

[3] A. Gathmann and M. Kerber, A Riemann-Roch theorem in tropical geometry. Math. Z. in press, 
preprint: math. AG/0612129 (2006). 

[4] R. Hirota and S. Tsujimoto and T. Imai, Difference scheme of soliton equations, in Future directions 
of nonhnear dynamics in physical and biological systems, eds. P.L. Christiansen, J.C. Eilbeck, and 
R.D. Parmentier (Plenum Press, New York, 1993). 

R. Hirota and S. Tsujimoto, Conserved quantities of a class of nonlinear difference-difference equa- 
tions, J. Phys. Soc. Japan, 64, no. 9, 3125-3127 (1995). 
[5] T. Kimijima and T. Tokihiro, Initial-value problem of the discrete periodic Toda equations and its 

ultradiscretization. Inverse Problems 18 1705-1732 (2002). 
[6] A. Kuniba, T. Takagi and A. Takenouchi, Bethe ansatz and inverse scattering transform in a periodic 

box-ball system, Nucl. Phys. B 747, no. 3, 354-397 (2006). 
[7] J. Mada, M. Idzumi and T. Tokihiro On the initial value problem of a periodic box-ball system. J. 

Phys. A 39 L617-L623 (2006). 
[8] G. Mikhalkin, Enumerative tropical algebraic geometry in R2, 'math. AG/0312530' 
[9] G. Mikhalkin and I. Zharkov, Tropical curves, their Jacobians and theta functions, math. AG / 0612267] 
[10] P. van Moerbeke and D. Mumford, The spectrum of difference operators and algebraic curves. Acta 

Math., 143, 93-154 (1979). 
[11] A. Nagai, T. Tokohiro and J. Satsuma, Ultra- discrete Toda molecule equation, Phys. Lett. A, 244, 

383-388 (1998). 



[12] D. Speyer and B. Sturmfels, Tropical Mathematics, |math.CO/0408099 

[13] D. Takahashi and J. Matsukidaira, On discrete soliton equations related to cellular automata, Phys. 
Lett. A, 209, 184-188 (1995). 

21 



[14] D. Takahashi and J. Satsuma, A soliton cellular automaton, J. Phys. Soc. Japan, 59, 3514-3519 (1990). 
[15] T. Tokihiro, D. Takahashi, J. Matsukidaira and J. Satsuma, From Soliton Equations to Integrable 

Cellular Automata through a Limiting Procedure, Phys. Rev. Lett., 76, 3247-3250 (1996). 
[16] F. Yura and T. Tokihiro, On a periodic soliton cellular automaton, J. Phys. A; Math. Gen. 35 , 3787- 

3801 (2002). 

Department of Physics, Graduate School of Science, The University of Tokyo, 7-3-1 
hongo, bunkyo-ku, tokyo 113-0033, japan 

CREST, JST, 4-1-8 Honcho Kawaguchi, Saitama 332-0012, .Japan 
E-mail address: reiiyOspin.phys . s .u-tokyo . ac . jp 

Faculty of Marine Technology, Tokyo University of Marine Science and Technology, 
2-1-6 Etchu-jima, Koto-Ku, Tokyo, 135-8533, Japan 
E-mail address: takeiiawa@kaiyodai.ac.jp 



22 



